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Motivated by recent progress on the understanding of the Eguchi-Kawai (EK) volume equivalence 
and growing interest in conformal window, we simultaneously use the large-A^^ volume reduction 
CN and Chiral Random Matrix Theory (;^RMT) to study the chiral symmetry breaking of four dimen- 

o 

sional SU(A^c) gauge theory with adjoint fermions in the large Nc limit. Although some cares are 
required because the ;^RMT limit and 't Hooft limit are not compatible in general, we show that 
the breakdown of the chiral symmetry can be detected in large-A^^ gauge theories. As a first step, 
we mainly focus on the quenched approximation to establish the methodology. We first confirm 
^ that heavy adjoint fermions, introduced as the center symmetry preserver, work as expected and 

thanks to them the volume reduction holds. Using massless overlap fermion as a probe, we then 
^ calculate the low-lying Dirac spectrum for fermion in the adjoint representation to compare to 

^ that of ;^RMT, and find that chiral symmetry is indeed broken in the quenched theory. 



The XXXth International Symposium on Lattice Field Theory - Lattice2012 
July 24-29, 2012 
Cairns, Australia 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos.sissa.it/ 



Large Nc volume reduction and chiral random matrix theory 



Jong-Wan Lee 



1. Introduction 

In modern theoretical particle physics, numerical simulations of quantum field theory on the 
lattice play an important role in non-perturbative studies of strongly coupled gauge theories, e.g. 
quantum chromodynamics (QCD), from the first principle. However, lattice calculations often 
encounter practical challenges such as large finite- volume effects and sign problems. For instance, 
lattice studies of the infrared behavior of conformal theories are extremely difficult due to finite 
size of the lattice. These challenges may be overcome by constructing lattice formulations for 
large-A^c gauge theories (for a review, see [1]) or for effective field theories [2], where various 
nice properties hold. In particular, the large-A^c volume equivalence [3] has recently received large 
attention because of the theoretical understandings of the center symmetry stabilization and the 
success of the numerical tests [4, 5, 6, 7, 8].^ 

One of the motivation of our work stems from the growing interest in conformal window 
and the walking technicolor model (WTM), where conventional lattice QCD techniques can be 
easily extended to QCD-like theories. A candidate of the minimal WTM is two-flavor SU (2) 
gauge theory with fermions in the adjoint representation [10]. Bearing this in mind, we consider 
numerical calculations of SU (Nc) gauge theories with Uf = 2 adjoint fermions in the large Nc limit. 
We first clarify the difference between the ;^RMT limit and 't Hooft large-A/^ limit in large-A/^ 
gauge theories. Then, we establish the methodology - the way to demonstrate the spontaneous 
chiral symmetry breaking (S;^ SB) by applying the ;^RMT techniques and the large-A^c volume 
equivalence to large-A/^ gauge theories (for earlier work along this direction, see [11, 12]). As a 
concrete example, we study the heavy adjoint two flavor QCD up io Nc = 16 on a 2"^ lattice, which 
approximates the quenched QCD, and find that the chiral symmetry is indeed broken. 

2. The 't Hooft versus ;f RMT limit in large-A/^ gauge theory 

In large-A/c gauge theories, the 't Hooft limit is the large-A^c limit in which the 't Hooft coupling 
A = g^Nc and the number of quark flavor^/ are fixed [13]. The space-time volume V and the quark 
mass m are also fixed. In this limit, the theory is dramatically simplified, i.e. the l/Nc expansion 
has a natural topological structure and a certain class of Feynmann diagrams, so called ''planar 
diagrams'', provide dominant contributions in the perturbation theory. For example, the vacuum 
expectation value of a properly normalized single trace operator S' can be expanded as, 

oo 

E CH,Bi^,V,m)N;^''-^, (2.1) 

where h and B are the number of handles and boundaries of diagrams. The connected correlation 
functions of more than one operators have the same structure. 

In the 't Hooft limit, various nice properties hold; in particular, vacuum expectation val- 
ues of products of gauge-invariant operators are factorized up to 0(1 /Nc) corrections, (G\G2) = 
((^i)(<^2) + 0(1/ Nc). An important consequence of the large-A^c factorization is the so-called 
Eguchi-Kawai (EK) volume equivalence [3]; provided center synmietry is not broken sponta- 
neously, the Wilson loop amplitudes in an arbitrarily small volume agree with those of the usual 
large-volume lattice gauge theory. The large-A/^ volume equivalence can also be understood as 
an example of orbifold equivalence [5] and the general statement is as follows. Starting with an 
original (parent) theory, we obtain a new (daughter) theory by performing a projection under some 

^Unbroken center symmetry is the most crucial condition for the large-A^c volume equivalence. The original proposal 
for pure U {Nc) gauge theory by Eguchi-Kawai failed due to spontaneous breaking of center symmetry [9]. 
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discrete subgroup of the global symmetry of the parent theory. If the discrete symmetry does not 
break down spontaneously, the correlation functions of invariant (neutral) sectors of operators in 
both theories are equal up to a trivial rescaling factor. 

Provided the chiral symmetry is spontaneously broken, the dynamics of QCD at below Aqcd 
can be described by the low-energy effective theory where the degrees of freedom are Goldston 
bosons such as pions instead of quarks and gluons. If we go further inside the e-regime, where 
the one-fourth of space-time volume is much smaller than the pion Compton length, the relevant 
degrees of freedom are zero modes of pions and the density of Dirac eigenvalues near zero is related 
to the chiral condensate Z via Banks-Casher relation [14]. The distribution of the low-lying Dirac 
eigenvalues can be calculated using the chiral random matrix theory (;f RMT) in the microscopic 
limit, N ^ oo while having mN fixed, where the matrix size N is identified by V in the standard 
gauge theory (for a review, see [15]). In large-A^c gauge theories, the number of color charges A^c 
are also relevant degrees of freedom and the corresponding limit may be obtained by taking the 
large-A^c with a fixed mN ^ mVN^, where the constant a can depend on the representation 
of fermions. Let us call this limit as the ";^RMT limit" which is different from the 't Hooft limit: 
the fermion mass scales with Nc and hence Ch^B in Eq. 2.1 also has a nontrivial A^c-dependence 
implying that the 't Hooft power counting rules and the large-A^c equivalences are not valid. In 
the next section, however, we argue that the EK equivalence and the ;f RMT combined with lattice 
simulations can be used to detect the S;^ SB in QCD-like theories with some cares. 

3. Strategy of detecting the spontaneous chiral symmetry breaking 

In this section we establish the way to use the ;f RMT in large-A/^ gauge theory for detection 
of the S;^SB. First let us recall how one can confirm the chiral symmetry breakdown in the or- 
dinary 5(7(3) QCD. The criterion for the SxSB is the nonzero chiral condensate in the standard 
thermodynamic limit, 

{xfry/) 7^ in the massless limit m^O after taking the large-volume limit y ^ oo. 

On the other hand, for QCD in the e-regime, the S;^SB is recognized by the agreement of the 
low-lying Dirac spectrum with the prediction from ;f RMT. 

In the case of large-A^c gauge theory, the logic is as follows. First, we calculate the low- 
lying Dirac spectrum of the large-A/^ gauge theory in a small box (e.g. on a or 2"^ lattice) and 
compare the spectrum with the prediction from ;^RMT. Here the probe mass ruprobe must scale 
as fUprobe ^ ^/^c large N^ otherwise such that the ;fRMT limit is realized. Because of this 
scaling, the EK equivalence, requiring Nc with m fixed, can not be applied directly. However, 
following the same logic as the QCD case except that V is now replaced with A^^, the agreement of 
the spectrum with the prediction from ;f RMT still establish the S;^SB of the large-A/^ gauge theory 
in a small box. Then the EK equivalence leads to the S;^SB of the large-A^c gauge theory at large 
volume. The power a is determined so that rUprobe and the Dirac eigenvalues near zero have the 
same Nc -dependence, and we will see a = 1 in our setup. The schematic diagram of the detection 
of the S;^SB in the large Nc limit is shown in Fig. 1. 

4. Numerical simulation of the Eguchi-Kawai model 

Numerical simulations are performed with the Wilson gauge action and Wilson-Dirac fermions 
in the adjoint representation on a 2^ lattice. 




(4.1) 
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Figure 1: Detection of the spontaneously breakdown of chiral symmetry in large-A/^ gauge theories using 
the ;^RMT and EK equivalence. 



and 

2 / 4 , 



J=l n \ Ai 



(4.2) 



where b and K represent the inverse of the 't Hooft coupling, Z? = 1 / g^Nc, and the hopping paremeter 
which is related to bare quark mass by K* = {2ma + 8)~^ respectively. The Plaquettes P^v are built 
of link variables in the fundamental representation as usual. For the fermionic action, the link 
variables in the adjoint representation are defined by 

Kb = \^AT§UTIu\ (4.3) 

where Tp are SU {Nc) generators in the fundamental representation. This action is invariant under 
SU {Nc) local gauge transformation as well as Z^^ global center transformation: 

C/n,M f2nC/n,Mf^n+M' ^neSU{Nc), and C/„,^ e^^'^'^/'^^Un^^, % G Z^V,- (4.4) 

Throughout our studies, we use periodic boundary conditions for all lattice directions in both link 
variables and fermion fields. 

Our lattice simulations consist of two parts: 1) quenched calculations (k* = or equivalently 
ma is infinite), as a nontrivial check of our numerical code by confirming the breaking of center 
symmetry at weak coupling, 2) simulations for two adjoint fermions whose mass is of order I /a 
(k = 0.09) ^, where low-lying Dirac eigenvalues are calculated by using a massless overlap-Dirac 
fermion as a probe. We performed simulations at = 0.5 for up to A^c = 16 which is relatively 
smaller than that used for simulations of a single-site EK model [6, 7, 16]. As we will see below, 
however, we could obtain good large Nc limits since we have additional suppression of the finite 
volume effects thanks to the larger volume y = 2^. For A/^ = 8 and K* = 0, we also performed 
simulations at = 0.3 and 0.4 corresponding to the strong and intermediate couplings, respectively. 
For all lattice simulations, we used the Hybrid Monte Carlo (HMC) algorithm; 200 trajectories are 
used for the initial thermalization, typically 500 configurations are generated for each ensemble, 
and every two adjacent configurations are separated by 10 trajectories. 



^These heavy adjoint fermions serves as a ''center symmetry preserver'', where our chose of K" = 0.09 and b = 0.5 
is inside of the conjectured center symmetric region in the K — b plane [16], but do not play any role in other parts of 
dynamics. Thus the case 2) can be also considered as a quenched calculation. 
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Figure 2: (Left) Scatter plots of the Polyakov lines, which are overaged over all 2^ sites orthogonal to the 
line dirction /i = 1, for A^^ = 8 and /? = 0.5. The red and blue dots are about 150 data points obtained from 
ensembles for /c = 0.09 and k* = 0, respectively. (Right) Ensemble averaged values of magnitudes of the 
Polyakov lines versus for the EK model coupled to two heavy adjoint fermioins. The red solid line 
represent the fits of the 5 < A/^ < 16 data to the function cq + ci /A/^. 



4.1 center symmetry 

In order to use the large-A^c volume equivalence, it is essential to confirm that the center sym- 
metry is unbroken. Some evidences of which center symmetry is unbroken are as follows: (1) the 
Polyakov line scatters radially in the vicinity of origin in the complex plane, (2) the magnitude 
of approaches zero as Nc increases where the predicted asypmtotic scaling behavior is 1 /Nc, (3) 
the average plaquette value measured from the reduced model agrees with that measured from the 
large- volume lattice gauge theory. The Polyakov line along /i -direction in a 7^ lattice is defined by 

P^(x) = -^TrC/^(x)C/^(x + ^^). (4.5) 

From numerical simulations for quenched EK model with various values of Nc and b, we first 
reproduced the well known results for pure Yang-Mills theory: the center synmietry is unbroken at 
strong coupling, but broken at weak coupling. In particular, we show scatter plots of the P^'s for 
A/c = 8 and b = 0.5 in Fig. 2, where the red and blue dots represent the data points for K = 0.09 
(EK model coupled to two heavy adjoint fermions) and /c = (quenched EK model), respectively. 
The two colors have very different behaviours: the red develops a cluster around the origin, which 
is consistent with the senario of unbroken center symmetry, while the blue is localized at some 
of the elements of the center of SU{Nc), which implies that the center symmetry is broken. At 
strong coupling, i.e. b = 0.3, the P^'s radially scatter around the origin in both simulations. We 
also calculated ensemble averages of for the case of /c = 0.09 and b = 0.5. As shwon in 
Fig. 2, they are linear in 1/A^c and extrapolate to zero at Nc = oo. Finally, we caculated the average 
plaquette values and performed extrapolations to Nc = oo. The values we obtained are 0.72053(71) 
for K = 0.09 and 0.72733(12) for K* = 0. Compared with the large-volume lattice simulation [6] 
and the single-site EK model for quenched QCD [7], the former is in good agreement while the 
latter is systematically larger. In summary, the center symmetry is intact in the theory with two 
heavy adjoint fermions and thus the EK volume equivalence holds; the euqivalent large- volume 
lattice theory approximates the quenched large-A/^ QCD. 

4.2 Dirac spectrum and comparison to ;f RMT 

Our strategy of detecting the S;^SB is to compare microcopic properties of the EK model with 
those of ;^RMT, where the ;^RMT limit without losing generality might be achieved by simply 
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Figure 3: (Left) Low-lying Dirac spectrum for A^^ = 16, = 0.5 and /c = 0.09, where z is the rescaled 
Dirac eigenvalue obtained by matching the ensemble average of the lowest eigenvalue with the expectation 
value of the lowest eigenvalue in ;^RMT for the ChGSE. Colored bars represent the histograms of individual 
eigenvalues, while the empty bars represent the accumulated histogram. (Right) Dirac eigenvalue spacings 
mutiplied by A/^ for A^^ = 8, 12, 16, where 5Xk = ^k — ^k-\ and Aq = 0. 



taking ruprobe = and Nc oo. For this purpose, we calculate the low-lying spectrum of the 
hermitian overlap-Dirac operator & for a massless fermion in the adjoint representation.^ The 
operator ^ is defined by 

^ = \r5 [1 + 75 sgn[i/^(-mo)]] , (4.6) 

where YsHw{—mo) is the Wilson-Dirac operator and the Wilson mass is set to be negative —mo < 0. 
Throughout this work, ;^RMT predictions are restricted to quenched (Nf = 0) theory with zero- 
topological charge. The adjoint QCD belongs to the universal class of the Chiral Gaussian Sym- 
pletic Ensemble (ChGSE). However, we also consider two other universal classes, Chiral Gaus- 
sian Orthogonal Ensemble (ChGOE) and Chiral Gaussian Unitary Ensemble (ChGUE), in order to 
make the comparison manifest. As shown in the left figure of Fig. 3, the distribution of low-lying 
Dirac eigenvalues is in good agreement with that of ChGSE, impling that the chiral symmetry is 
spontaneously broken. 

It might be interesting to find how rUprobe scales with A^c to extend our study to non-zero mass 
of the probe fermion or to determine the size of the effective volume in the EK model. In the 
right figure of Fig. 3, we plot the spacings of low-lying Dirac eigenvalues (ensemble averaged) 
multiplied by Nc for Nc = 8, 12, 16. From the figure, we see a nice agreement between the spacings 
at up to ^ = 2 for A^c = 8 and at up to ^ = 3 for Nc = 12. This agreement implies that the Dirac 
eigenvalues near zero, which are expected to agree with ;f RMT preeiction, scale as Nc and thus 
a = 1. We also found that the spectral density p(z) is nearly zero at between the {Nc — l)th and 
A/cth eigenvalues. 

A possible explanation of our founding for the A^c-scaling of Dirac eigenvalues may rely on 
the perturbative analysis in the background of diagonal Wilson lines [7]. In a compact space and 
at weak coupling, one cannot gauge away zero-momentum modes and thus the low-lying Dirac 
spectrum might be determined by the zero modes. For fermion in the adjoint representation, the 
number of zero modes of the Wilson lines is {Nc — 1) while the total degrees of freedom is N^. As 
a result, the low-lying Dirac eigenvalues should scale as N^ which is consistent with our numerical 
results. Interestingly, our counting of zero modes also agrees with the position at which p(z) ^ 0. 

^The Dirac eigenvalues for adjoint fermions are appearing as conjugate pairs with 2-fold degeneracy and we take 
the positive distinct eigenvalues for our numerical results. 
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5. Conclusion and outlook 

In this proceeding, we discussed about the 't Hooft limit and ;^RMT limit in large-A^c gauge 
theories: in general two limits are not compatible and the large-A^c equivalences do not hold in 
;^RMT limit. In spite of the difference of these two limits, the S;^SB can be detected by taking 
an indirect path from xRMY limit of EK model to the large-volume lattice theory in the 't Hooft 
limit. As a numerical demonstration, we performed lattice simulations of SU {Nc) gauge theory 
with two heavy adjoint fermions on a 2^ lattice. After confirming that the volume equivalence is 
valid (unbroken center symmetry), we found that the distributions of low-lying Dirac eigenvalues 
are in good agreement with the ;^RMT prediction and thus the chiral symmetry of the quenched 
QCD is spontaneously broken. In the near future, we hope to determine whether the large-A/^ 
dynamical two-flavor adjoint QCD goes through the S;^SB or not, with an application to walking 
Technicolor theories in mind. 
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